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DIAGNOSTIC APPARATUS AND METHODS 
FOR A CORIOLIS FLOW METER 

Background of the Invention 

5 

1. Field of the Invention 

The present invention relates to diagnostic apparatus and methods for a Coriolis flow meter. 

2. Statement of the Problem 

10 It is known to use Coriolis mass flow meters to measure mass flow and other 

information of materials flowing through a pipeline as disclosed in U.S. Patent No. 
4,491,025 issued to J.E. Smith, et al. of January 1, 1985 and Re. 31,450 to J.E. Smith of 
February 11, 1982. These flow meters have one or more flow tubes of different 
configurations. Each conduit configuration may be viewed as having a set of natural 
15 vibration modes including, for example, simple bending, torsional, radial and coupled 
modes. In a typical Coriolis mass flow measurement application, a conduit configuration is 
excited in one or more vibration modes as a material flows through the conduit, and motion 
of the conduit is measured at points spaced along the conduit. 

The vibrational modes of the material filled systems are defined in part by the 
20 combined mass of the flow tubes and the material within the flow tubes. Material flows 
into the flow meter from a connected pipeline on the inlet side of the flow meter. The 
material is then directed through the flow tube or flow tubes and exits the flow meter to a 
pipeline connected on the outlet side. 

A driver applies a force to the flow tube. The force causes the flow tube to oscillate. 
25 When there is no material flowing through the flow meter, all points along a flow tube 
oscillate with an identical phase. As a material begins to flow through the flow tube, 
Coriolis forces cause each point along the flow tube to have a different phase with respect to 
other points along the flow tube. The phase on the inlet side of the flow tube lags the 
driver, while the phase on the outlet side leads the driver. Sensors are placed at different 
30 points on the flow tube to produce sinusoidal signals representative of the motion of the 
flow tube at the different points. A phase difference of the signals received from the sensors 
is calculated in units of time. 

The phase difference between the sensor signals is proportional to the mass flow rate 
of the material flowing through the flow tube or flow tubes. The mass flow rate of the 
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material is determined by multiplying the phase difference by a flow calibration factor. 
Prior to installation of the flow meter into a pipeline, the flow calibration factor is 
determined by a calibration process. In the calibration process, a known fluid is passed 
through the flow tube at a given flow rate and the proportion between the phase difference 
and the flow rate is calculated. 

One advantage of a Coriolis flow meter is that there are no moving components in 
the vibrating flow tube. The flow rate is determined by multiplying the phase difference 
between two points on the flow tube and the flow calibration factor. The phase difference is 
calculated from sinusoidal signals received from the sensors indicating the oscillation of 
two points on the flow tube. The flow calibration factor is proportional to the material and 
cross sectional properties of the flow rube. Therefore, the measurement of the phase 
difference and the flow calibration factor are not affected by wear of moving components in 
the flow meter. 

However, it is a problem that material properties, cross sectional properties and the 
stiffness of a flow tube can change during use of the Coriolis flow meter. The changes in the 
material properties, cross sectional properties and stiffness of the flow tube are caused by 
erosion, corrosion, and coating of the flow tube by material flowing through the flow tube, 
changing pipeline mountings and temperature. One example of the change in cross-sectional 
properties of the flow tube is the change in the moment of inertia caused by corrosion of the 
flow tube. A second example of a change in the material and cross-sectional properties of 
the flow tube is an increase of the mass of the flow tube and a decrease in cross-sectional 
areas caused by coating of the flow tube by materials flowing through the tube. A change in 
the material properties, cross sectional properties and stiffness of the flow tube can change 
the flow calibration factor of the flow meter. If the flow calibration factor of the flow meter 
changes, flow rates that are calculated using the original flow calibration factor are 
inaccurate. Therefore, there is a need in the art for a system that detects a possible change in 
the material properties, cross sectional properties and/or stiffness of a flow tube indicating 
that the mass flow rates measured by the Coriolis flow meter may be inaccurate. 
Summary of the Solution 

The above and other problems are solved and an advance in the art is achieved 
through the provision of a system for validating the flow calibration factor of a Coriolis 
flow meter through the deterrnination of flexural stiffness and other parameters. The 
described validation system can be used to solve either single degree of freedom models or 
multiple degrees of freedom models. 
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Some examples of the validation system capable of solving single degree of freedom 
models include measuring the tube stiffness with an auxiliary displacement sensor. The 
auxiliary displacement sensor is placed on a flow tube of the flow meter. A known force is 
applied to the flow tube and the resulting deflection is measured by the auxiliary 
5 displacement sensor. The resulting tube stiffness can be use to predict changes in the sensor 
calibration factor. 

Some examples of the validation system capable of solving single degree of freedom 
models include estimating tube stiffness from a drive loop frequency response. The drive 
loop response function contains a stiffness term which can be correlated to the flexural 
10 stiffness of the flow tube. The resulting flexural stiffness can be used to predict changes in 
the sensor calibration factor. 

Some examples of the validation system capable of solving multiple degree of 
freedom models include estimating tube stiffness from model parameters. A response model 
representing the flow meter is converted to a modal model. The modal model is than 

i 

15 converted to a physical model. The physical model contains parameters of interest 
including tube stiffness and damping that can be used to determine a change in the 
calibration factor. 

In some examples, the data acquired must be normalized. Normalization is needed 
to cancel out changes in the properties of the system that change with time. Correcting 
20 resonant frequency due to changes in tube temperature is one example of the need for 
normalization. 

Some examples of the validation system include methods and appartus for 

correction of the flow calibration factor using linear or non-linear equations, trending or 

multi-fluid calibration. 
25 Description of the Drawings 

FIG. 1 illustrates a Coriolis flow meter in an example of the invention; 

FIG. 2 illustrates a validation system in an example of the invention; 

FIG. 3 illustrates a validation system in an example of the invention; 

FIG. 4 illustrates a validation system in an example of the invention; 
30 FIG. 5 illustrates a validation system in an example of the invention; 

FIG. 6 illustrates a validation system in an example of the invention; 

FIG. 7 illustrates a validation system in an example of the invention. 
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Detailed Description of the Invention 

FIGS. 1-7 and the following description depict specific examples to teach those 
skilled in the art how to make and use the best mode of the invention. For the purpose of 
teaching inventive principles, some conventional aspects have been simplified or omitted. 
Those skilled in the art will appreciate variations from these examples that fall within the 
scope of the invention. Those skilled in the art will appreciate that the features described 
below can be combined in various ways to form multiple variations of the invention. As a 
result, the invention is not limited to the specific examples described below, but only by the 
claims and their equivalents. 
Example 1 

Factors that affect flexural stifBiess also affect Coriolis flow meter sensitivity (flow 
calibration factor). Flexural stiffiiess is the static spring rate derived from flexing the flow 
tube with a known force pattern and measuring the flow tube displacement. Any force 
pattern could be used to measure flexural stiffiiess, as long as it is invariant. As an example, 
the flexural stiffiiess for a clamped beam is as follows: 

K Flex = F = 192EI (1) 
8 L 3 

where: 

F - Force (N); 

E - Young's Modulus (N/m 2 ); 

I - Moment of Inertia (m 4 ); 

L - Length (m); 

Kfiex - flexural stiffiiess of flow tube. 

For a Coriolis flow meter, if flexural stiffiiess changes, then so must calibration 
factor. Flexural stiffiiess of a Coriolis flow meter is defined as: 

Kflex = CpC G Cs[EI] (2) 

where: 

C P - effect offeree pattern on flexural stiffiiess; 

Cg - effect of unflexed tube bend geometry on flexural stiffiiess; 

Cs - effect of unflexed tube stress on flexural stiffiiess. 

For a straight tube Coriolis flow meter with no pre-stress the following expressions 
show the dependence of calibration factor on EI: 

m=c[^] AT- (3) 
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So the flow calibration factor (FCF) for the straight tube is: 

fcf= 4f ] w 

where C is a constant determined by mode shape and pick-off locations. 

FIG. 1 illustrates a Coriolis flow meter 5 capable of directly estimating flexural 
stif&ess by applying a known force to the flow tubes and measuring the resulting 
deflection. Coriolis flow meter 5 comprising a flow meter assembly 10 and meter 
electronics 20. Meter electronics 20 is connected to meter assembly 10 via leads 100 to 
provide density, mass flow rate, volume flow rate, totalized mass flow and other 
information over path 26. 

Flow meter assembly 10 includes a pair of flanges 101 and 101', manifold 102 and 
flow tubes 103A and 103B. Connected to flow tubes 103 A and 103 B are driver 104, pick- 
off sensors 105 and 105' and auxiliary position sensor 107. Brace bars 106 and 106' serve to 
define the axes W and W about which each flow tube 103A and 103B oscillates. 

When flow meter assembly 10 is inserted into a pipeline system (not shown) which 
carries the material being measured, material enters flow meter assembly 10 through flange 
101, passes through manifold 102 where the material is directed to enter flow tubes 103 A 
and 103B, flows through flow tubes 103 A and 103B and back into manifold 102 where it 
exits meter assembly 10 through flange 101'. 

Flow tubes 103 A and 103B are selected and appropriately mounted to manifold 102 
so as to have substantially the same mass distribution, moment of inertia, and elastic 
modulus about bending axes W~W and W'-W respectively. The flow tubes extend 
outwardly from the manifold in an essentially parallel fashion. 

Flow tubes 103A-B are driven by driver 104 in opposite directions about their 
respective bending axes W and W f and at what is termed the first out of phase bending 
mode of the flow meter. Driver 104 may comprise one of many well known arrangements, 
such as a magnet mounted to flow tube 103 A and an opposing coil mounted to flow tube 
103B. An alternating current is passed through the opposing coil to cause both tubes to 
oscillate. A suitable drive signal is applied by meter electronics 20, via lead 110 to driver 
104. The force pattern generated by driver 104 is fixed over the life of the sensor (i.e. C P is 
constant). Driver 104 generates a constant force and auxiliary position sensor 107 
determines the tube displacement. As shown in FIG. 1, the auxiliary sensor would be 
positioned near the pick-offs 105 or 105' however other suitable locations could be utilized. 
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Meter electronics 20 receives the right and left velocity signals appearing on leads 
111 and 111 1 , respectively and the tube displacement over lead 112. Meter electronics 20 
produces the drive signal on lead 110 which causes driver 104 to oscillate flow tubes 103 A 
and 103B. The present invention as described herein can produce multiple drive signals 
5 from multiple drivers. Meter electronics 20 processes left and right velocity signals and a 
tube displacement signal to compute mass flow rate and provide the validation system of the 
present invention. Path 26 provides an input and an output means that allows meter 
electronics 20 to interface with an operator. 

FIG. 2 shows a process 200 for validating the flow calibration factor of Coriolis 

10 flow meter 5 using calculated changes in its flexural stiffness. Process 200 begins in step 
210 by applying a known force to flow tubes 103 A&B using driver 104. In step 220, 
auxiliary position sensor 107 determines the tube displacement resulting from the force 
applied by driver 104. The flexural stif&iess, Knew, is measured in step 230 using the force 
and tube displacement from steps 210 and 220. A previous flexural stiffness, Kow, is 

15 retrieved in step 240. This previous flexural stiffness is determined either at the factory or 
after the meter is installed. Step 250 compares K* C w with Koid to determine if there has been 
a change in the flexural stiffness of the meter. If no change has been detected, step 260 
signals a "no error'* condition and repeats the process. If a change has been detected, step 
270 corrects the flow calibration factor (FCF) based on the stif&iess change. 

20 Example 2 

Flow tube flexural stiffiiess can also be determined by estimating points on a tube 
frequency response function (FRF) at given frequencies. These points are then used to fit a 
single degree of freedom model to the data and finally determine the DC (e.g. zero crossing) 
point on the FRF. 

25 It is possible to write a general expression for any cross (not collocated) mobility 

frequency response function as a linear combination of SDOF responses: 



or equivalently modeled by: 

30 
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where v is the pickoff velocity, /is the driver force and H(s) is the FRF of interest. The 
drive loop frequency response can be approximated well in this manner. In fact, at 
frequencies near the drive frequency, only one term in the summation is required for good 
approximation of the frequency response. If more accuracy is required in parameter 
estimation, then more than one mode can be retained in the model. 

As an alternative to including more modes, the higher frequency modes could be 
masked from the drive loop frequency response by selecting the location of the driver and 
pickoff coils. The second bend mode (twist mode) can be masked from the drive loop 
frequency response because the driver is located on the node of the second mode. If the 
pick-ofFs were located on the nodes of the third mode, then the third mode is masked from 
the drive loop frequency response as well. With this driver/pickoff pattern, the drive loop 
frequency response is very close to a SDOF system. 

Modeling the drive loop as an SDOF system results in the following: 

^H(s)= ' (7) 
f m,5 + c x s + k x 



or equivalently modeled by: 

— = H(s) = i—ttt T ( 8 > 

The drive loop frequency response can be converted from mobility (velocity) to 
receptance (displacement) by integrating (multiplying by 1/s ): 

x = H ^ - 1 (9) 
/ s rr^s 2 + Cj5+ k x 

again equivalently modeled by: 

x _ H(s) = 4 (jo) 

F s s 2 + 2£ x & x s+ co 2 

The stiffiiess term k\ is well correlated with flexural stiffiiess, and thus is a good predictor 
of changes in flow cal factor caused by changes in flexural stiffiiess. 

To estimate the stiffiiess parameter k\ a single or multiple frequency estimation can 
be used. FIG. 3 shows a process 300 for validating a flow calibration factor using single 
frequency estimation. Single frequency estimation 300 begins in step 310 by determining a 
receptance transfer function (9). In step 320, the inverse of the receptance transfer function 
at ia> is determined to get the inverse receptance frequency response: 

-1!-^ — L = m, (i(D) 2 + c x ic>+ k x (1 1) 

H(l(D) 1 
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Next in step 330, we let ©—►(): 

lim tttt : = limwtfioj) 2 + limc 1 ia>+ limfc, (12) 



lim J ^ = lim/:, = (13) 

<y->o H{lCO) Q>^>Q 1 



So for co— >0: 



(14) 

1 H{ia>) 



If mi is known or is negligible, then a larger frequency closer to co\ can be used to estimate 
10 k\ from: 

k x - r^rr - qto (15) 
If ci and mi are both known, then any frequency can be used to estimate k\ from: 

* ! " H^) 1ThiiQ>) " ( } 

In step 340 we calculate M from expressions (14), (15), or (16) with co in radians/second and 
1 5 correct the flow calibration factor based on the calculated k\ in step 350. 

FIG. 4 shows a process 400 for validating a flow calibration factor using multiple 
frequency estimation. Multiple frequency estimation 400 begins in step 410 by identifying 
constants mi, ci, ku cou and Ai using any time domain or frequency domain system 
identification method. In step 420, a curve fitting procedure is used to fit a rational 
20 continuous time transfer function model to the complex frequency response vector H at the 
set of frequencies in vector W (in radians/second). The number and location (in frequency) 
of the FRF data points does affect the quality of the fit. A good fit is achieved using as few 
as 2 frequency response data points. The derived model is of the form: 

b(J)s* + b(2)s< N >-» + -+KN b + \) 



The above model is incompatible with the mobility (velocity) form of the drive loop 
frequency response, because there is no way to force b(Nb + 1) = 0. In step 430, we get 
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around this limitation by converting the driver pickoff mobility (velocity) frequency 
response data to the receptance (displacement) form. The measured mobility frequency 
response data if must be multiplied by \l{i(6). The measured mobility drive loop frequency 
response H should be from drive coil current (proportional to force) to pickoff voltage 
(proportional to velocity). 

Converting the mobility data to receptance data yields H(s) in the form: 

^ )= a(iy + a(2> + a(3) (18) 
where a(l) = L In step 440, the modal parameters of interest are extracted from the transfer 
function model as follows: 

4 = b(i) 

a>, = R3) (19) 

The physical parameters can then be calculated in step 450 using the following equations: 
77*! = 1 / A } 

c x = 2£ x G>JA, (20) 

Once the physical parameters are determined, changes in the flow calibration factor 
as well as other parameters, including changes in the mass and length of the flow tube, are 
determined and corrected in step 460. In addition to correction, warnings and meter shut 
down can be initiated based on the predetermined criteria. 
Frequency Response Function Determination Using Multi-sine Excitation 

Accurately estimating points on the FRF can be difficult in real-time. Windowing 
and spectral line estimation often have estimation bias which has a direct impact on the 
ability to predict FCF changes. In order to more efficiently determine estimating points on 
the FRF and eliminate the need for windowing, a multi-sine input is used. 

The multi-sine input method applies an efficient spectral line estimation routine 
using a minimum crest factor multi-sine input excitation signal. The minimum crest factor 
signal is used to minimize peak drive power and to limit large output signals. Large output 
signals often have non-linear components that cause measurement errors. 

FIG. 5 shows a process 500 for calculating spectral lines used in validating a flow 
calibration factor using multiple frequency estimation. Process 500 begins with step 510 by 
determining measurement frequencies of interest. The frequencies of interest are 
determined from the following: 

9 
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f j+1 =fj + AfO=l:N sin ) (21) 
Next, a multi-sine excitation signal is defined in step 520 as: 

AM Nsin-1 

* = £ £ sin(2^// ; + ^) (22) 

I«0 JmQ 

5 where, 

fj = frequencies for multi-sine 
(|>j = phase angle for multi-sine where -n < <|>i < n 
ti = discrete time point 
N = number of data points in multi-sine 
1 0 N S in = number of sine waves in multi-sine 

In order to maximize the input energy while minimizing the peak value, a Crest Factor (CF) 
is minimized with respect to <|>j in step 530. Crest Factor (CF) is defined as max / rms, 

CF(s)- , maX( * (0) (23) 

V J (AM 

x fZx(if/N 

15 

In order to make efficient use of the multi-sine, total measurement time, step 540, and 
total number of averages, step 550, must be determined. The total measurement time can be 
calculated using the following equation: 

tp = 1/ Af (24) 
20 Equation (24) yields a measurement time equal to the period of the multi-sine, resulting in a 
measurement that is minimally affected by parameter changes. The total number of 
averages required is normally chosen so as to reduce the standard deviation of the estimate 
due to noise. The signal-to-noise ratio and its direct relationship to the estimated 
confidence interval, normally determine the total number of averages necessary for a 
25 required accuracy. 

Step 560 applies the multi-sine, jc, to the input of the system and measures the 
resulting output, y. Step 570 uses x and y to determine the transfer function, stiffness and 
calibration factor using the techniques described herein. 
Example 3 

30 A changing flow calibration factor can also be detected and corrected by utilizing a 

multiple degree of freedom model. In a multiple degree of freedom model, a change in the 
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flow calibration factor is a function of model parameters (i.e. mass, stiffness, mode shape, 
complex pole, etc.), fluid properties (density, pressure, etc.) and environmental conditions 
(temperature, vibration etc). 

In order to determine parameters that affect the flow calibration factor, a response 
5 model must first be determined. Response models are valuable in that they are direct 
measurements of the actual structure. The structure's response from a known excitation is 
usually in the form of an impulse response function (IRF) or frequency response function 
(FRF). After the response model is determined, the response model is converted into a 
modal model, and then into a spatial model. 
10 For example, in a modal model, the flow calibration factor is a function of the mode 

shape and modal frequencies as evidenced by the following expressions: 
L%MFCF- f{ Pjluid ,K)= C fPjluid + C K K 
A%MFCF=g(0,A)=C o 0 + C X X 

where 

A%MFCF - percent change in mass flow calibration factor 
15 p fl Uid - fluid density 

K - stiffness of the sensor. 
O r . mode shapes 

- complex pole of the sensor's r th mode. 
By periodically measuring modal parameters, fluid properties and other environmental 
20 parameters, a change in the flow calibration factor can be detected and corrected. 

FIG. 6 shows a process 600 for validating the flow calibration factor of a Coriolis 
meter using a multi-degree of freedom model. Process 600 begins in step 610 by creating a 
response model of a structure. Step 620 converts the response model from step 610 into a 
modal model. The modal model is then converted into a spatial model in step 630. In step 
25 640, physical parameters, such as flexural stiffness, can be determined from the spatial 
model. 

Step 650 compares the current physical parameter data with previously determined 
physical parameter data, to establish if a change has occurred. If no change has occurred, 
measurement can continue to take place at a predefined interval at the request of a user. If a 
30 change has occurred, an error signal is generated in step 660. After the error signal in 
generated, step 670 corrects the calibration factor. Correction can either be done 
automatically or only at the request of the user. 

11 
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Correction of the flow calibration factor can be accomplished using various methods 
including coefficient estimation, multi-fluid calibration and trending. Coefficient 
estimation uses equations, like equation 25 above, with an estimation of C f and C K . For 
instance, if the material properties of a meter are affected by temperature, values for Cf and 
Ck can be estimated by finite element analysis and determination of a nominal material 
property change with temperature. 

Multi-fluid calibration involves the calibration of a meter using at least two fluids at 
two different stiffnesses. Coefficients Cf and Ck can be determined by using the following 
least squares fit: 



MFCF, 

M 

MFCF, 



Pfluid x 
M 

Pjluidfi 



M 



(26) 




P fluids 
M 



K 



N 



MFCF, 

M 

MFCF N 
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Trending techniques can also be used to correct the flow calibration factor. Initial 
model parameters, fluid properties, and temperature data is collected in the transmitter and 
the coefficients for a polynomial, n-dimensional "operating surface 11 f(pjj uid ,K,T,o;P 9 ...) 

would be determined. Because there is generally a linear relationship between flow 
calibration factor, fluid properties, model parameters, and temperature, a change would be 
identified when newer data showed coefficients changing with respect to the initial 
coefficients. 

As an example, assume the fluid density has no effect on flow calibration factor, or 
Cf is zero. The measured stiffness will change with temperature due to changes in the flow 
tube's material properties. Furthermore the measured stiffness will change due to 
temperature affecting a transducer's material properties and sensitivity. Therefore the rate of 
change of stiffness due to temperature is unknown. Nevertheless, the stiffness and 
temperature can be measured at a specific time, i.e. time of installation. The data can be 
used to establish the relationship between stiffness and temperature. Most material 
properties change linearly with temperature over a wide range as shown by: 

K(T) = /(!> C^A 7-+ (27) 

where 

Ck = rate of change of stiffness due to temperature, 
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C trans i = rate of change of the transducers sensitivity change, and 

Ctrans2 = transducer's offset 
Initial data taken at the time of installation or at the factory over a range of flow tube 
temperatures allows you to calculate Ck, Ctnmsi and Ctrans2. 
5 The stiffness measured at any temperature will fall in a trend line. When enough 

stiffness measurements are made over a reasonable temperature range, the slope, or the rate 
of change of the stiffness with respect to temperature can be estimated using a least squares 
fit. This is the common change of the stiffness with respect to temperature. When the 
stiffness and temperature are measured and the data point does not fall on the curve, then a 
10 change in the sensor has occurred and the MFCF has also changed. 

The fact that the stiffness estimate is changing with temperature indicates that the 
flow tube material and transducer material properties are in slight error. A correction can be 
made in the temperature correction coefficients to make the measured stiffness vs. 
temperature slope zero. The correction may be made by an optimization process to 
15 minimize the error. Subsequent to this correction, a stiffness change will indicate the change 
in the MFCF. 

The concept of a multiple degree of freedom model can be applied to other 
parameters, including flow calibration factor changes due to fluid density, pressure, external 
vibration, etc. 
20 Normalization 

The described solutions to a multiple-degree of freedom model are ideal when the 
structure does not change over time. However, in some applications a structure changes 
over time. In the case of a flow sensor, the density of a process fluid conveyed by the 
sensor may change over time resulting in a change in mass of the sensor. Furthermore, 
25 ambient temperature and process fluid temperature may change over time, altering the 
material properties of the sensor. In order to collect meaningful data at varying conditions, 
a normalization process needs be implemented. 

For the case of a flow sensor, there are three main types of normalization. The first 
is to normalize that data with respect to a resonant frequency. Various process fluid 
30 properties, such as fluid density, flow, and pressure, cause the resonant frequency to shift. 
The second is to normalize that data with respect to a reference temperature. The third is to 
normalize the data with respect to the displacement, or some other response variable. 

For example, a typical Coriolis flow meter operates at a specific resonance for a 

specific material density. At one fluid density, the resonant frequency is at 152 Hz. If the 
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fluid density increases, the resonant frequency shifts down to 121 Hz. The position of the 
resonant frequency relative to the test frequencies causes significant errors in the estimates 
of the model parameters. Therefore, the data needs to be normalized in order to produce 
meaningful results. 

5 FIG, 7 shows a process 700 for normalizing data in order to correct the flow 

calibration factor of a Coriolis meter using a multi-degree of freedom model. Process 700 
begins in step 710 by receiving data from the structure. Once the data is received it is 
normalized with respect to the resonant frequency. This allows the model parameters to be 
calculated from data measured for any fluid density. Tracking the resonant frequency allows 

10 the test frequency to be set as a scalar of the resonant frequency. 

Next in step 730, we normalize the data with respect to a reference temperature. This 
corrects errors resulting from data that is measured at two different temperatures. For 
example, two important material properties change with temperature causing a discrepancy 
in the measured data. The modulus of elasticity of the tube material changes by 

15 approximately -4.2%/100°C. Since most Coriolis flow meters use voice coil transducer to 
measure flow and drive, changes in temperature will change the strength of the magnetic 
field and change the length of the wire in the coil. For a typical AINiCo magnet, the 
magnetic flux density changes by -2.1%/100°C for both the driver and pickoff. 

Also, a Coriolis flow meter's modulus will change with temperature. Coriolis flow 

20 meters typically utilize a temperature sensor in order to report the temperature of the fluid in 
the sensor. Using this temperature measurement, the FRF measurements can be corrected 
back to some reference temperature. 

In step 740, the data is normalized with respect to displacement or some other a 
response variable. This is done by integrating from mobility (velocity/force) to receptance 

25 (displacement/force). The following is an example of the equation used to normalize the 
measured data: 



Once the data has been normalized, step 750 corrects the flow calibration factor using the 
30 methods described herein. 
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